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ABSTRACT 



o 

H 

GRB041219A is one of the longest and brightest gamma-ray bursts (GRBs) ever 
^ i observed. It was discovered by the INTEGRAL sateUite, and thanks to a precursor 

happening about 300 s before the bulk of the burst, ground based telescopes were able 
to catch the rarely-observed prompt emission in the optical and in the near infrared 
bands. 

Here we present the detailed analysis of its prompt 7-ray emission, as observed 
with IBIS on board INTEGRAL, and of the available X-ray afterglow data collected by 
\^ • XRT on board Swift. We then present the late-time multi-band near infrared imaging 

\ data, collected at the TNG, and the CFHT, that allowed us to identify the host galaxy 

• . of the GRB as an under-luminous, irregular galaxy of 5 x \Q^Mq at best fit redshift 

ofz = 0.3lt°Jt 



> 



o 



We model the broad-band prompt optical to 7-ray emission of GRB041219A 
within the internal shock model. We were able to reproduce the spectra and light 
curve invoking the synchrotron emission of relativistic electrons accelerated by a series 
of propagating shock waves inside a relativistic outflow. On the other hand, it is less 
easy to simultaneously reproduce the temporal and spectral properties of the infrared 
data. 

Key words: gamma-rays: observations - gamma-rays: bursts - galaxies: photometry 



1 INTRODUCTION 

Gamma-ray bursts (GRBs) were discovered in the late 1960s 
(Klebesadel et al. 1973), and 40 years later they are still not 
fully understood (see e.g. Gehrels et al. 2009, for a recent re- 
view). They are short lived transients (ms to hundreds of sec- 
onds) of soft gamma ray radiation that appear at random di- 
rections on the whole sky. A major breakthrough in GRB sci- 
ence has been achieved in the late 1990s, when thanks to the 
rapid localizations performed by the Italian-Dutch satellite 
^eppoSAX the first transient counterparts, the so-called af- 
terglows, have been discovered at X-ray (Costa et al. 1997), 
optical (van Paradijs et al. 1997), and radio wavelengths 
(Frail et al. 1997). These and the subsequent measurements, 
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thanks to the HETE-II, INTEGRAL, and Swift satellites 
have made it possible to firmly prove the cosmological na- 
ture of these sources, being their redshifts, z, distributed in 
the range (0.1-^8.2). Thanks to the distance measurements, 
the energy and luminosity of GRBs could be determined 
making them the most powerful explosions in the Universe, 
since their isotropically emitted energy iJiso spans from 10*° 
to 10*^ erg (Amati 2007). However, this huge amount of en- 
ergy request can be reduced in the hypothesis that GRBs are 
coUimated sources (Rhoads 1997), and indeed the detection 
of some achromatic breaks (e.g. Zeh et al. 2006) in the light 
curves of GRB afterglows further supports this interpreta- 
tion. From the measure of the break times the jet opening 
angles, Sjet, can be derived (Sari et al. 1999a), and taking 
this into account the energy reservoir and its spread can 
be reduced, clustering around 10*^ erg (Frail et al. 2001; 
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Bloom ot al. 2003). In several cases late bumps in the af- 
terglow light curves have been reported (e.g. Bloom et al. 
1999). These bumps have been interpreted as an emerging 
supernova component. In a handful of cases the GRB/SN as- 
sociation has been spectroscopically confirmed (see Woosley 
& Bloom 2006, for a review). These arguments and the en- 
ergetic budgets point towards the collapse of a massive (> 
30 Mq) star as being at the origin of GRBs, at least for the 
events lasting longer than ~2 s (Woosley et al. 1993). 

GRB 041219A was detected and localized in real time 
by the INTEGRAL Burst Alert System (IBAS; Mereghetti 
et al. 2003). The alert was issued when the burst was still 
on going, and this allowed some robotic telescopes to detect 
a prompt optical and infrared flash (Vestrand et al. 2005; 
Blake et al. 2005) whose position was consistent with the 
IBAS one (Gotz et al. 2004). 

Only a handful of bursts have been observed in the vis- 
ible range during the prompt phase, and, as often in GRB 
physics, it is impossible to identify a single common be- 
haviour (see Yost et al. 2007a,b). GRB 990123 was bright 
and showed no clear link between the three measurements 
made during the prompt phase and the gamma-ray data 
(Akorlof et al. 1999; Briggs et al. 1999). The "naked eye 
burst" GRB 080319B was even brighter, and may show some 
correlation with a delay of a few seconds between the op- 
tical and the high energy emission (Racusin et al. 2008). 
GRB 050820A (Vestrand et al. 2006) and GRB 041219A 
were much weaker (both in apparent and absolute magni- 
tudes) and their optical fluxes were apparently correlated 
to the 7-rays. In these two latter cases it seems reasonable 
to believe that the prompt emission can be explained by 
the same mechanism at high and low energy (Genet et al. 
2007b). Conversely, the two optically bright bursts proba- 
bly require different physical conditions or/and a different 
radiative process in the regions emitting the low and high 
energy photons respectively (Sari & Piran 1999b; Meszaros 
& Rees 1998; Zou et al. 2009; Hascoet et al. 2010). 

GRB041219A is the longest and brightest GRB in 
the IBAS sample (Vianello et al. 2009). Indeed the very 
good quality of the gamma-ray data allowed Gotz et al. 
(2009) to measure the time-variable high linear polariza- 
tion of this event in the 200-800 keV energy band, partly 
confirming previous but less significant results (Kalemci 
et al. 2007; McGlynn et al. 2007). Here we present the 
detailed spectral analysis of the prompt 7-ray emission of 
GRB 041219A, as measured by IBIS, and SPI-ACS on board 
INTEGRAL (Winkler et al. 2003), of its X-ray afterglow, as 
measured by Swift /XKT, and the identification of its host 
galaxy, as observed by the TNG, and the CFHT, as well as 
the modelling of the broad band, optical to 7-ray, prompt 
emission of the GRB. 



2 OBSERVATIONS AND DATA REDUCTION 

2.1 INTEGRAL observations 

IBAS triggered at 01:42:17 UTC (from now on To) on De- 
cember 19*'' 2004 (IBAS Alert # 2073), on a precursor of 
GRB 041219A indicating a 38.2 sigma source at R.A. = 
24>" 26^ Dec. = -^62° 50' 06" (I = -M19.86°; b = 
-1-0.13°) with a 2.5' uncertainty . We selected the corre- 
sponding INTEGRAL pointing (SCWID: 026600780010), 



and processed the data using the Off-line Scientific Analysis 
(OSA) software provided by the INTEGRAL Science Data 
Centre (ISDC, Courvoisier et al. 2003) v8.0. 

We analysed IBIS {INTEGRAL coded mask imager; 
Ubertini et al. 2003) data. Our analysis is based on data 
taken with ISGRI (Lebrun et al. 2003), the IBIS low en- 
ergy detector array, which is made of CdTe crystals, and is 
working in the 13 keV-1 MeV energy range, and on PICsIT 
(Di Cocco et al. 2003) a pixellated Csl detector, working 
in the 150 kcV 10 MeV energy range. In the ISGRI 20-100 
keV image the only detected source is GRB 041219A, which 
is localized at R.A. = OO'' 24"" 27^ Dec. = -|-62° 50' 21" 
with a 19" uncertainty, consistent with the IBAS and the 
infrared (Blake et al. 2005) positions. We used this position 
to extract ISGRI spectra in 62 logarithmically spaced bins 
between 13 keV and 1 MeV. ISGRI light curves have been 
obtained by binning (At=l s) the ISGRI events in three 
different energy bands (20 40; 40-100; 100-300 kcV). The 
background has been subtracted by fitting a constant value 
to these light curves using the data of the same pointing 
before the bursts itself (i.e. before 01:41:00 UTC). 

The source is also clearly detected in the PICsIT single 
events histograms (252 336 keV band), but the positional 
uncertainty in this case is larger and we do not report the 
PICsIT position here. PICsIT histograms are not suited for 
spectral extraction in this case, since they are integrated 
over the entire duration of the jjointing (~1800 s). For spec- 
tral and light curves extraction we used PICsIT spectral 
timing data. They have no positional information, since they 
represent the total count rate of the PICsIT camera in 4 en- 
ergy bands (158-208; 208-260; 260-364; 364-676 keV). We 
used them to extract light curves, rebinned at 4 s intervals. 
These light curves have been background subtracted by eval- 
uating the background before and after (next pointing) the 
GRB. The same data have then been used to extract PICsIT 
spectra by summing the count rates relative to the GRB in 
4 bands and rebiiining the PICsIT response matrix in order 
to cope with available energy bands. Due to the fact that the 
GRB was in the fully coded field of view, no further correc- 
tions were necessary, and the standard ancillary response file 
could be used. For the spectral analysis, background spec- 
tra have been derived by using the pointing before and the 
pointing after the GRB. 

For timing purposes only, we also used the data of the 
Anti-Coincidence System (ACS) of the INTEGRAL spec- 
trometer SPI (Vedrenne et al. 2003). ACS data consist of 
the total count rate (50 ms time resolution) above ~100 
keV measured by the 91 Bismuth Germanate (BGO) scin- 
tillator crystals that surround SPI. The crystals are used 
as the anti-coincidence system of the spectrometer, but are 
also capable of detecting high-energy transient events such 
as bright GRB and SGR bursts (von KienUn et al. 2003). 
We computed the ACS light curve with a bin size of 5 s and 
estimated the background by fitting a constant value to the 
data of the same pointing excluding the burst itself. 

2.2 Swift/XRT observations 

At the epoch of GRB 041219A, the automated slewing of 
the Swift narrow- field instruments was not enabled yet: for 
this reason the X-Ray Telescope (XRT) did not begin col- 
lecting data until about 5 hours after the IBIS detection. 



GRB 041 21 9 A and its host galaxy 3 



Swift /yiRT observed the afterglow for roughly 0.6 ks in win- 
dowed timing (WT) mode and for about 1 ks in low rate 
photodiode (LR) mode (sec Hill et al. 2004, for more details 
on the XRT data-taking modes). Given the low count rate, 
here we only considered WT data. 

The data were processed with standard procedures and 
screening criteria using the xrtpipeline (version 0.12.3 un- 
der the Heasoft package 6.6). We extracted the source 
events from a 40 x 40 pixels (1 pixel corresponds to about 
2.36 arc sec) box along the image strip; a region of the same 
shape and size, positioned well outside the point-spread 
function of the source, was used to estimate the background. 
The ancillary response file was generated with xrtmkarf, 
and it accounts for vignetting and point-spread function cor- 
rections. 

2.3 TNG/NICS observations 

The first attempt to identify the host galaxy of GRB 
041219A, was made through an H-band observation of the 
field of the GRB with the NICS camera mounted on the 3.5 
m TNG^ telescope at La Palma. 

Observations of the field of GRB041219A were per- 
formed on Oct 15"* 2007, i.e. a few years after the GRB. 
The observations were secured under good observing condi- 
tions with about 1 " seeing in the H band. In total 15 frames 
with 30 s exposure and NDIT=3 were collected. Data reduc- 
tion was carried out following standard procedures with the 
Eclipse package (Devillard 2001). The median of the input 
frames was derived to obtain a sky frame which was subse- 
quently subtracted from the input frames. Then sub-pixel 
registration was applied, and a final scientific frame weis ob- 
tained averaging all input frames. An astrometric solution 
was derived by a large number of 2MASS'^ objects and the 
final absolute accuracy is about 0.5" both in R.A. and Dec. 
Aperture photometry was computed with tools provided by 
the GAIA package^ and the absolute calibration was derived 
by comparison with a suitable number of not-saturated, iso- 
lated 2MASS stars in the field. 

2.4 CFHT/WIRCam observations 

Further observations of the field of GRB 041219A were ob- 
tained during the night of October 6*'' 2009 (program ID; 
09BF08) with the WIRCam instrument at the 3.6 m CFHT 
in Mauna Kea under good sky conditions (seeing 0.5-0.8"), 
providing a complete set of photometric images. The region 
of GRB 041219A was observed in the Y, J, H, and filters 
for about 30 minutes per filter. 

In order to perform sky subtraction, a first stack of the 
images has been produced using a median of the input im- 
ages. Saturation masks have been generated looking for sat- 
urated pixels in the individual images, in order to discard 
saturated objects from the catalogues. A global astromet- 
ric solution has been found for the whole dataset using the 
2MASS catalogue as an absolute reference. A single binary 
mask, masking all the sources detected with sextractor has 

^ http://www.tng.iac.es/ 

^ http://www.ipac.caltech.edu/2mass/ 

^ http://star-www.dur.acuk/~pdraper/gaia/gaia.html 



then been produced, and the sky value for each pixel could 
be determined. On the sky subtracted images relative pho- 
tometry and astrometry could then be performed. The inter- 
nal dispersion (difference between same objects coordinates 
in different images) of the astrometric solution for high S /N 
objects is small, ^0.06", and the reference dispersion (dif- 
ference between the object coordinates in our images and 
the 2MASS catalog ones) of the astrometric solution for the 
high S/N objects is around 0.13" (similar to the internal 
accuracy of the 2MASS astrometry) . 

The stability of the sky emission was not perfect during 
this run. As a consequence, a number of images show strong 
large scale gradients in the background, even after the sky- 
subtraction, because of the short time scale fluctuations of 
sky emission lines. The following procedure was applied to 
remove this pattern before the combination of the images: 
we ran sextractor on the individual images with a mesh 
size of 256 pixels and using a very aggressive masking of 
the objects (produced during the sky-subtraction step). The 
output background image was then subtracted to the input 
image. 

One stack is finally produced for each of the four fil- 
ters {Y , J, H, and Ks). The absolute photometry is then 
checked: while for the J, H, and Ks bands the 2MASS cata- 
log has been used as a photometric reference, for the Y band 
zero points have been adopted from CFHT observations of 
standard stars'*. 

3 ANALYSIS AND RESULTS 
3.1 INTEGRAL 

3.1.1 Timing analysis 

GRB041219A, being the brightest and longest burst de- 
tected by INTEGRAL, is the only one for which good PIC- 
sIT data arc available. Therefore we included them in our 
analysis. In this section we present the analysis of this GRB 
using both IBIS detector planes and the Anti-Coincidence 
System (ACS) of the INTEGRAL spectrometer. 

The light curve of GRB041219A, as measured by the 
different high-energy instruments is shown in Fig. 1. The 
GRB consists of two main peaks, which are preceded by 
two precursors, the second of which is detected only in the 
ISGRI data. The ISGRI data are affected by telemetry sat- 
uration for what concerns the first precursor, and the two 
main peaks, hence these light curves are not representative 
of the true ISGRI recorded flux. PICsIT and ACS data, also 
shown in Fig. 1, are complementary to the ISGRI ones in 
the sense that they extend the observable energy range to 
higher values, but, in view of the globally smaller number 
of registered counts, they are not telemetry limited on what 
concerns the two main peaks. 

3.1.2 Spectral analysis 

ISGRI spectra have been derived for the four different sec- 
tions of the GRB identified as first and second precursor, and 

* http: / /www.cfht. hawaii.edu/Instruments/Imaging/ 
WIRCam/IiwiVersionlDoc.html#Part6 
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Figure 1. Background-subtracted light curve of GRB 041219A as seen by ISGRI, PICsIT, and the Anti-Coincidence System of SPI in 
different energy bands. The error bars have been omitted for clarity. 



first and second peak, as defined in Table 1. PICsIT spectra 
were obtained by integrating the spectral timing data over 
the same time intervals considered for the ISGRI analysis. 
Only for the first and second peak PICsIT data have enough 
statistics to be included in the spectral fits. 

The data were grouped with a minimum of 20 counts 
per energy bin and the spectra were analysed with the xspec 
fitting package version 11.3. A multiplicative factor has been 
added to the joint fits in order to take into account the in- 
struments inter-calibration. The two datasets never needed 
to be adjusted by more than 30%. For fitting purposes a sim- 
ple power law model has been considered as a first approach, 
then an exponential high-energy cut-off has been adopted to 
better describe the data, and finally a "Band" GRB model 
has been used if necessary. The results of the fits are re- 
ported in Table 1. Due to the high fiux, 3% of systematics 
has been added to the spectra during the fit. 

The first precursor has a very hard spectrum with re- 
spect to the rest of the GRB, as is seen also from the light 
curves inspection (see Fig. 1). It can be well fit by a single 
power law model, but a slight improvement can be achieved 
if the spectrum is fitted with a cut-off power law model, 
see Table 1. Even if the cut-off energy, Ec, is not well con- 
strained for what concerns its maximal value, a solid lower 
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Photon Index 

Figure 2. Confidence contours (68%, 90% and 99% c.l.) of the 
cut-off energy versus the photon index F for the ISGRI fit of 
the first precursor. 



limit can be derived, as shown by the contour plot (photon 
index vs. cut-off energy) of Fig. 2. The derived (90%) lower 
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Figure 3. Simultaneous ISGRI/PICsIT spectral fit of the first 
peak of GRB 041219A. In the upper panel the data points (crosses 
for ISGRI and triangles for PICsIT) are compared to the best fit 
model, while the lower panel represents the residuals with respect 
to the fit. 

limit for the cut-off energy {E^) is 180 keV , corresponding 
to an Speak energy^ larger than 130 keV, and a best fit value 
of 370 keV, making this part the hardest of the whole GRB. 

The second precursor shows a soft spectrum, and it can 
be described by a simple power law model without the need 
of any cut-off. Other models like a black body or the quasi- 
thermal model (see below) can be ruled out (x^/d.o.f. — 
343.5/46, and x^/d.o.f. = 57.3/42, respectively). It is the 
softest part of the whole GRB, as it is not detected above 
100 keV. 

The first main peak is very bright, and in this case the 
single power law fit can be excluded (xVd-o.f. = 116.63/59), 
while the cut-off power law model and the Band model give 
an equally good representation, with x^/d.o.f. = 72.29/58 
and /A.o.i. — 72.11/57. But having the Band model an 
extra degree of freedom, and being the parameter /3 not well 
constrained, in Table 1 we report only the fits done using 
the cut-off power law model. The simultaneous ISGRI and 
PICsIT best fit model and data are shown in Fig. 3. 

The second peak is fainter and softer than the first one. 
As for the first peak, it can equally well be fit by a Band 
model or by a cut-off power law, and the best fit of the 
joint ISGRI/PICsIT data is reported in Table 1. The simple 
power law model, on the other hand, is excluded (x^/d.o.f. 
= 103.97/54). The derived -Epcak in this case is significantly 
smaller than the one derived for the first peak, as can be 
seen by comparing the two values in Table 1. 

The average ISGRI/PICsIT spectrum of the two main 
peaks is also well fit by a cut-off power law, and its param- 
eters are also reported in Table 1. Its fluence in the 20-200 
keV energy band is 2.5x10"* erg cm~^. The bolometric (1 
keV-10 MeV) fluence of the whole GRB is ~ 5x10"* erg 
cm^^. The derived spectral parameters can be used to de- 
termine the pseudo — z (Atteia 2003; Pelangeon et al. 2006) 
for GRB041219A, which resuhs to be 0.19±0.06. Using this 
distance estimator the isotropic equivalent energy released 
by this GRB would be of the order of 2x10^^ erg. 

We note that, except for the first precursor, the spectral 

^ -Bpeak is given by (2-r)xEc in this case. 
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Figure 4. Time resolved spectral analysis of the IBIS/ISGRI 
data. The derived photon index values are plotted against time 
(lower panel), and the 20-300 keV light curve is plotted in the 
upper panel, as a reference. 

parameters derived here are compatible within the errors 
with the ones derived by McBreen et al. (2006), obtained 
using the germanium spectrometer, SPI, on board INTE- 
GRAL. The differences for the first precursor, that in our 
analysis is softer, may result from the much better statistics 
we have at low energies in the current analysis. Furthermore 
using the quasi-thermal model, composed by the sum of a 
black body and a power law, for modelling the rest of the 
GRB, as proposed by the same authors, we cannot confirm 
their findings in terms of black body temperatures, obtain- 
ing fcT = 16 ± 6 keV for the first peak and fcT = 19 ± 5 keV 
for the second, while they obtain much higher values. 

Thanks to the fact that GRB 041219A is very bright, we 
have been able to perform a detailed time-resolved analysis. 
We extracted contiguous spectra integrated over 20 s during 
the main GRB emission (i.e. excluding the precursors). Due 
to the lower signal-to-noise ratio, the spectra were just fitted 
with simple power-law models. The results are presented in 
Figure 4. As can be seen, a classical hard-to-soft evolution is 
observed, and we note that the photon index F evolves from 
about 1.5 to about 2.5, indicating that while -Epcak decreases 
during the burst (as indicated by the values in Table 1) it 
literally crosses the IBIS/ISGRI energy domain, since the 
values of F measured at the end of the burst are likely to 
be the asymptotic value of the /3 parameter of the Band 
function. 

3.2 X-ray afterglow 

The observed net count rate during the WT mode observa- 
tions was as low as 0.60 ± 0.03 counts s~^ between 0.3 and 
10 keV, so no pile-up correction was necessary. The data 
were grouped with a minimum of 20 counts per energy bin. 
Given the paucity of counts, we fit a simple model to the 
data: a power law corrected for absorption. We obtained the 
following best-fitting (x^ ~ 1-0 for 14 degrees of freedom) 
parameters: absorbing column (1.5 ± 0.3) x lO^'^ cm~^ and 
photon index F = 1.9 ± 0.3. The observed (not corrected 
for absorption) 0.3-10 keV flux was ~ 4.ll°;^ x 10"" erg 
cm"'^ s"^ at To + 5 hours, while the absorption corrected 
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Table 1. INTEGRAL spectral parameters. Errors are given at 90% c.l. 



Name 



Tstart 
UTC 



Tstop 
UTC 



keV 



-^peak 

keV 



d.o.f PICsIT 



First Precursor 


01:42:17 


01:43:12 


1.41±0.04 






60.52 


49 


no 




01:42:17 


01:43:12 


1.26±0.13 


>180° 


>134'' 


57.67 


48 


no 


Second Precursor 


01:43:27 


01:44:12 


2.16±0.06 






31.41 


44 


no 


First Peak 
Second Peak 
1=' peak+ 2"'' peak 


01:46:22 
01:48:12 
01:46:02 


01:47:40 
01:48:52 
01:50:42 


1.42±0.05 
1.71±0.06 
1.53±0.01 


347+ J.f 
334+if 
353+1^^ 


201+80 


68.37 
78.03 
77.59 


58 
52 
62 


yes 
yes 
yes 



The best fit values are 370 keV and 274 keV for E^ and i?peaki respectively. 
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Figure 5. CFHT/WIRCAM ii's-band image of the region around 
GRB041219A. The circle represent the position of the GRB in- 
frared flash, as derived by Blake et al. (2005). The black labels 
correspond to the three objects discussed in the text. 



All uncertainties 



one was ~ S.Ql^^o ^ 10~^^ erg cm~^ s 
are quoted at la confidence level. We note that the absorb- 
ing column density we derive is larger than the value re- 
ported by Kalberla et al. (2005) (~8xl0^^ cm"^), but the 
two values are compatible at a ~2 a level, and the resolu- 
tion of the LAB HI survey, being of the order of 0.6°, is too 
coarse for clumpy low-latitude Galactic regions. But the dif- 
ference may also be accounted for by the absorption within 
the GRB host galaxy, being this excess value typical of GRB 
afterglows (e.g. Campana et al. 2010). 



3.3 Near Infrared Photometry 

We have observed the field of GRB041219A, obtaining first 
H-hand photometry with the NICS camera at the Italian 
Telescopic Nazionale Galileo, as well as later photometry in 
the YJHKa bands with WIRCam mounted at the Canada- 
France Hawaii Telescope. One object is clearly detected in 
our images with a position compatible with the one of the 
prompt infrared flash (Blake et al. 2005). 

Data corresponding to this object (labeled "Object 3", 
see Fig. 5) are hsted in Table 2. We only list the CFHT 
results, whose higher quality supersedes the TNG data. We 
must remark that the position of this source is within the 
Galactic plane (|ti| < 1°), and thus it is expected to be heav- 
ily reddened. 



Parameter 


Value 


R.A. (J2000) 


00'' 24™ 27.6' 


Dec. (J2000) 


-1-62° 50' 33. 5" 


Y 


22.16 ±0.35 


J 


20.81 ±0.20 


H 


19.73 ±0.17 




18.86 ±0.12 



Table 2. Data corresponding to the GRB host, as obtained 
through our CFHT observations. All magnitudes are given in the 
Voga-based Johnson system, and are not corrected for Galactic 
extinction. 



3.3.1 Object Identification 

We have compared the YJHKs colours of Object 3 to sev- 
eral models of different galaxy types, in order to obtain some 
information about the probable redshift and spectral type. 
In Fig. 6 we present the results of such comparison. We 
have performed the same analysis for two nearby objects 
(labeled 1 and 2, and already noticed by Blake et al. 2005), 
but they were not compatible with typical colours of galaxies 
of any type. The objects have been de-reddened using the 
local value as given in Schlegel et al. (1998), which yields 
E{B - V) = 1.858 and Ay = 5.760 mags - the correction 
is also shown in the plots. Given the low Galactic latitude 
of the GRB and the coarse resolution of the Schegel maps 
(2°), we checked the absorption value against our X-ray af- 
terglow observations applying the relation between optical 
extinction and hydrogen column density derived by Giiver 
& Ozel (2009), A^H {cm~^) = (2.21 ±0.09) x 10^^ Ay (mag). 
Using the Ah value derived from our X-ray afterglow obser- 
vations, we obtain Av ~ 6.8 ± 1.6 mags, which is compatible 
within the errors with the value derived by Schlegel et al. 
(1998), that has hence been used as the reference value in 
this work. 

Figure 6 shows that the de-reddened colours of Object 3 
are compatible with moderate redshift solutions. We have in- 
deed used the photometric redshift code by Fernandez-Soto 
et al. (1999) to estimate the redshift and basic properties of 
the host galaxy, finding that the best fit solution corresponds 
to an irregular galaxy template at redshift z — O.SIIqj^. We 
must remark, however, that there is a high degeneracy of the 
models in that region of colour space, which renders possi- 
ble a large subset of galaxy models, spanning from early- 
type galaxies at low redshifts (z ~ 0.1 to 2 ~ 0.5) to late- 
type models at moderate redshifts {z ~ 0.2 to z 1.8). 
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galaxy, together with our best-fit spectrum, eorresponding to an 
' ' - irregular galaxy at redshift z = 0.31. We also include in the plot 
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Figure 6. Four-band YJHK colour-colour diagrams for the af- 
terglow host. In both panels we present the AB-magnitude colour 
tracks for different galaxy templates, ranging from elliptical to 
irregular as indicated in the legend therein. Each track begins at 
redshift z = 0, in the position marked with a filled circle, close 
to colour = 0. From there the redshift increases, passing through 
z = 1 (empty circle), z = 2 (cross), and reaching 2 = 3 (filled 
square). The observed and de-reddened colours of the GRB host 
are marked with triangles, surrounded by boxes which indicate 
the associated 1-sigma errors. We have added in quadrature an 
extra error in the de-reddened colours to include the uncertainty 
in the reddening value, assuming a 10% uncertainty in each term. 

Indeed by taking all the possible galaxy types into account 
the best fit redshift is 2 = 0.3, with a bimodal la confi- 
dence interval in the ranges 0.1-1.1 and 1.3-1.7 (secondary 
minimum). The uncertainties are estimated by evaluating 
the observed photometric uncertainties, that are included 
in the fitting procedure, and the systematic ones, induced 
by the inability of the spectral templates to include all the 
variability present in galaxy spectra. Fernandez-Soto et al. 
(2002), using BVRIJHK data, showed that the latter could 
be modelled with a normal error distribution with a vari- 
able sigma a{z)/{l -I- z) « 0.065. Given our similar set- 
up, but the diflerent wavelength coverage, we estimate that 
(j[z) / {1 + z) ~ 0.10 is an adequate assessment of the system- 
atic error. Using this estimate and convolving the two error 
terms as in Fernandez-Soto et al. (2002) we calculated the 
confidence interval mentioned above. 

However, there are other pieces of information that help 
in solving this degeneracy. In particular, an application of 
the Amati relation (Amati et al. 2002) clearly prefers the 
redshift range z < 0.5, and the pseudo-redshift technique 
(see 3.1.2) yields a value z — 0.19 ± 0.06^. In addition, re- 

^ Those two indicators are not independent, as the Amati relation 
is one of the factors taken into account in the calculation of the 
pseudo-redshift as described by Atteia (2003). 



garding the spectral properties, an early-type galaxy is, ac- 
cording to previous experience (e.g. Savaglio et al. 2009), 
less likely to be the host of a long GRB. 

Taking all this information into account, we have 
adopted the best-fit z — 0.31 as reference value for the red- 
shift, and an irregular template as the spectral choice. With 
these values, we show in Fig. 7 the intrinsic and observed 
spectrum of such an object. As can be seen, the photometry 
fits the observations perfectly, and we can derive an abso- 
lute rest-frame K-hand magnitude K = —21.0 AB, which 
corresponds approximately to a ~ 0.1 L* (Fontana et al. 
2006) galaxy, and, using a simple relation (e.g. Savaglio et 
al. 2009), to a mass of 5 x 10^ M©. The derived redshift 
yields a luminosity distance of 1.6 Gpc (for Ho = 71 km s~^ 
Mpc~^) and a total isotropic emitted energy Eieo ~ 10^'' 
erg (1 keV-10 MeV). However, taking the uncertainty on 
the photometric redshift into account the allowed intervals 
for the quantities cited above are [—16.5, —23.6] for the 
K AB mag nitude, [3 x 10^ 7 x 10^°] M© for the mass, and 
[2 X 10^\ 2 X 10^''] erg for Ei^o- 

3.4 Broad-band SED modelling of the prompt 
emission 

3.4-1 Constructing the broad-band SED 

Thanks to the rapid response to the IBAS alert of the RAP- 
TOR telescope, Vestrand et al. (2005) were able to mea- 
sure the optical emission of GRB 041219A during its prompt 
phase. There are three R-band data points and two upper 
limits that indicate that the optical emission is variable and 
likely correlated to the 7-ray one (see Fig. 3 of their paper). 
Blake et al. (2005) recorded the emission of GRB 041219A in 
three infrared bands, but just on the very final part of the 
prompt 7-ray emission. We use these visible and infrared 
data in conjunction to those at high energy presented in 
this work to build the broad-band SED of the prompt GRB 
emission at different times. 
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To do this we extracted four 7-ray spectra correspond- 
ing to the three time periods where the optical emission has 
been recorded, and to the one where the infrared data are 
available. In Table 3 we report the selected time intervals 
we used and the results of the fits. The infrared data have 
been de-reddened using the same parameters as in Section 
3.3.1. 

For the first three intervals a cut-off power law model 
was needed to fit the 7-ray data, while for the last one no 
cut-off could be measured, and as stated before, we may 
be measuring the l3 parameter of a Band function (see Sec. 
3.1.2). 



3.4-2 Modelling the prompt emission of GRB 04.1219A 

To model the SED of the prompt phase of GRB 041219A 
from the infrared to the 7-ray range we have considered the 
standard scenario of internal shocks (Rees & Meszaros 1994) 
assuming that the observed radiation comes from the syn- 
chrotron emission of shock accelerated electrons (Piran 1999; 
Bosnjak et al. 2009). In order to reproduce the complex 
time history of the burst we have supposed that the initial 
Lorentz factor distribution in the relativistic outfiow is made 
of a succession of elementary injections. The kinetic power 
is not the same for all injection phases but remains constant 
during one given episode (see top panel of Fig. 8) . The to- 
tal injected kinetic energy is Skin = 1-95 x 10®'* erg. The 
adopted micro-physical parameters, which are necessary to 
estimate typical values for the post-shock electron Lorentz 
factor and local magnetic field, are : ee = 0.33, es = 0.0033 
and C = 10~^ (fraction of electrons which are accelerated). 
The spectrum for each individual shock is computed using 
the standard synchrotron prescriptions of Sari et al. (1998). 

Each elementary injection episode develops its own sys- 
tem of internal shocks and a fraction €e of the dissipated 
energy is then radiated by the electrons resulting in a simu- 
lated gamma-ray light-curve represented in the middle panel 
of Fig. 8 (see Daigne & Mochkovitch 1998, for a detailed de- 
scription of the model). We do not try to get an exact fit of 
the observed light curve, but it can be seen in Fig. 8 that 
the different phases and corresponding levels of the prompt 
activity are correctly reproduced. 

The simulated spectra for the first three time inter- 
vals are also in reasonable agreement with the observa- 
tions (Fig. 8, bottom panel) while the infrared data dur- 
ing the final interval appear more difficult to reproduce (see 
Sect. 3.4.4 below). The adopted slope for the population of 
accelerated electrons is p = 2.8 in order to get the correct 
photon spectral index in the fourth interval, assuming that 
it indeed represents the j3 parameter of a Band function. 
In the other intervals we have no way to constrain 13. We 
have then kept the same value p = 2.8. The maximum of 
each spectrum (in yF^) is located at the synchrotron fre- 
quency Vrn while the evolving cooling frequency Vc remains 
located between and the optical domain. This allows us 
to reproduce the different optical levels in the first three 
time intervals. Also note that the self-absorption frequency 
fa always remains well below the optical band. 



3.4.3 Deceleration by the external medium 

The discussion presented in the last section only takes into 
account internal shocks and neglects the deceleration effect 
by the burst circumstellar medium. This is justified as long 
as the typical radius of internal shocks is smaller than the 
deceleration radius. In the case of GRB 041219A, which had 
a very long duration and varied over a wide range of time 
scales (from less than one second to several tens of seconds), 
this condition implies that the burst environment must have 
a low density : n < 0.1 cm~^ for a uniform medium and 
At < 0.01 for a stellar wind where A, is the wind param- 
eter A = M/4ttv normalized to 5 x 10^^ g cm~^ which is 
the value expected for a standard Wolf-Rayet wind having 
M = 10"® M0 yr-i and v = 1000 km s"* (Chevalier & 
Li 1999; Eldridge et al. 2006). These limits are quite re- 
strictive, and, if the density is larger, the reverse shock will 
propagate inside the jet before the end of the prompt phase, 
therefore affecting the dynamics. 

Hence we tried to model the temporal and spectral 
properties of GRB 041219A in the case where the decelera- 
tion effect of the external medium cannot be neglected. To 
follow the dynamical evolution of the ejecta we then adopt 
the method developed by Genet et al. (2007a). The reverse 
and internal shocks are now contemporaneous and their in- 
teraction modifies the observed emission; the location and 
intensity of the pulses in the 7-ray profile are changed, some 
are suppressed and others are created by late additional col- 
lisions in the decelerating flow. As a result, a new initial 
distribution of the Lorentz factor (Fig. 9, top panel) must 
be adopted to recover the main features of the temporal 
profile. The micro-physical parameters for the reverse shock 
are identical to those adopted for internal shocks. This as- 
sumption seems reasonable since the reverse shock, which is 
mildly relativistic and propagates within the ejecta, is very 
similar to internal shocks. 

Satisfactory results for the temporal and spectral prop- 
erties of GRB 041219A can then be obtained with a denser 
environment. Fig. 9 presents an example of a possible solu- 
tion for a wind like medium with A* = 0.1. 



3.4.4 Infrared emission 

The infrared flare detected during time interval 4 (see Table 
3) cannot be easily explained by internal shocks. Its soft 
spectral slope (photon index Fir ~ 1.7) differs from the 
gamma to optical average slope (F-^-opt ^ 1) and the limited 
data do not allow to assess its level of correlation with the 
high-energy emission. 

We therefore checked if it could be accounted for by a 
contribution from the forward shock. Tlion, the soft spec- 
tral slope can be reproduced but getting the correct time 
evolution appears very challenging. The late rise of the flare 
imposes too low values of the Lorentz factor and/or external 
density and moreover the observed variability of the infrared 
flux following flare maximum is not expected from the for- 
ward shock (see however Fan et al. 2005). 

The infrared flare may also come from a late collision in 
the flow but it is now not easy to get the right photon index. 
The observed Fir is close to 3/2 which would imply that Vc 
is below the infraxed domain whereas it was required to be 
above the optical at earlier times. Such a variation of Uc is 
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Table 3. Derived 7-ray spectral parameters, optical and infrared fluxes for the four intervals with simultaneous optical and infrared 
data. 



Selected Time 
Interval (UTC) 


Time since Tq 
s 


r 


keV 


-Epeak 

keV 


X^/d.o.f 


Flux density^ (filter) 
mJy 


01:45:43-01:46:56 
01:47:08-01:47:38 
01:47:51-01:49:04 


211 - 284 
296 - 326 
339 - 412 


1.3±0.1 
1.44±0.08 
1.79±0.07 


336+?i? 


175+if 
1881-8 

84lf3 


50.6/56 
50.6/56 
67.0/56 


2.9±0.8 {Rc) 
10.3±1.1 {Rc) 
3.8±0.8 (Rc) 



01:49 



16-01:51:54 



424 - 582 



2.41±0.03 



56.9/46 <1.1 (Rc); 1.7±0.1 (J); 2.0±0.1 (H); 2.6±0.1 (Ks) 



^Corrected for Galactic absorption (see text). Rc data points are taken from Vestrand et al. (2005), while the infrared data points are 

taken from Blake et al. (2005). 



not expected, except if the micro-physical parameters can 
vary by large factors during the propagation of the outflow. 

One could also try to solve the problem within a dif- 
ferent framework, where the infrared-optical and gamma- 
ray components result from different radiative mechanisms. 
For example in the Synchro-Self Compton scenario the low- 
energy photons correspond to synchrotron emission from 
relativistic electrons while the gamma-rays result from In- 
verse Compton scatterings of the synchrotron photons by 
the same electron population. This scenario has been con- 
sidered for the naked eye burst GRB 080319B (Racusin et al. 
2008; Zou et al. 2009; Kumar & Panaitescu 2008). However 
it has boon shown that in this case most of the dissipated 
energy is boosted in the GeV-TeV range by a second IC 
scattering (leading to an energy crisis, Piran ct al. (2009)). 
Moreover if this was the common scenario for GRBs, a sec- 
ond peak at high energy would then be expected for most 
GRBs, which has not been confirmed by Fermi observa- 
tions. Another possibility has been proposed by Zheng et al. 
(2006). In their model the gamma-rays penetrate an elec- 
tron cloud (produced in the supernova explosion prior to the 
GRB) and are partially converted into lower energy photons 
by successive Compton scatterings. 



3.4-5 Precursor activity 

As can be seen from Fig. 4, where the photon indices of 
all the spectra have been plotted versus time, there is a 
clear spectral evolution from hard to soft during the burst, 
superimposed on an intensity-hardness correlation, where 
the hardest spectra appear in correspondence with the peaks 
in the light curve. This "classical" evolution is however in 
contradiction with what is measured for the two precursors. 
If we take them into account, the hardest part of the GRB is 
represented by the first precursor and the softest one by the 
second precursor. Then the GRB becomes again very hard. 

Precursors have been subject to different definitions in 
the past. For instance Koshut et al. (1995) define precursors 
as events that are dimmer than the main GRB peak and for 
which there is a background interval at least as long as the 
main event. This does not strictly apply to GRB041219A. 
On the other hand, Lazzati (2005) identifies precursors as 
events that contain 0.1-1% of the counts with respect to 
the main event, and are separated by a period of quiescence 
without duration constraints. Our two precursors fit these 
requirements, since we find that they contain 1-2% and less 
than 1% of the total GRB flux, respectively. 

Precursors of thermal origin are expected in GRBs, due 



to the photo-spheric emission of the fireball when it becomes 
transparent (e.g. Daigne & Mochkovitch 2003), and may 
have been observed in the past (e.g. Murakami et al. 1991). 
However in our case neither the first nor the second precur- 
sor can be fitted with thermal spectra, and, in addition, the 
large delay between the precursors and the main emission is 
much longer than what is expected for thermal precursors. 

Different models have been invoked to explain non- 
thermal precursors. Waxman & Mcszaros (2003) and 
Ramirez-Ruiz et al. (2002) propose that they are associ- 
ated with the jet breakout of the stellar envelope, but the 
predicted delay between the precursor and the main GRB 
component is again too short to be applied to GRB 041219A 
(Wang & Mcszaros 2007). 

More generally, the ultra relativistic motion contracts 
so much the durations for the observer, that long delays are 
unreachable in all models where the same ejected material is 
supposed to radiate at different stages of its propagation to 
successively produce the precursor and the prompt emission. 
Therefore considering the similarity between the spectral 
properties of the precursors and the main emission - as also 
pointed out by Burlon et al. (2008, 2009) using BATSE 
and Sunft bursts -, we are left with a scenario of intermittent 
central engine activity, which remains to be explained in the 
context of the coUapsar model. 



4 CONCLUSIONS 

We presented a detailed spectral and temporal analysis of 
GRB041219A using IBIS and SPI-ACS on board INTE- 
GRAL. We have shown that the burst presents three spec- 
trally distinct phases (two precursors and a main event) that 
are separated by two long intervals of inactivity. We have 
analysed the available X-ray afterglow data of GRB 041219A 
constraining the column density in the direction of this burst 
located at low Galactic latitude. We performed a near in- 
frared observation campaign of the error region of the source, 
using TNG and CFHT, and identified the host galaxy of the 
GRB. Thanks to our multi-band photometry we modelled 
the host galaxy, and classified it as an under-luminous, likely 
irregular galaxy at rcdshift z = 0.3llg 26- Using this dis- 
tance estimate we were able to build the SED of the prompt 
emission of GRB041219A from the near infrared to the 7- 
rays. 

Our modelling of the light curves and broad-band spec- 
tra obtained for GRB 041219A shows that the prompt opti- 
cal and gammarray emission can be explained with a com- 
mon mechanism : synchrotron radiation from shock accel- 
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Figure 8. Model for the prompt emission of GRB 041219A. The 
top panel shows the adopted initial distribution of the Lorentz 
factor in the flow together with the value of the kinetic power for 
each injection episode as a function of the injection time tjnj . The 
middle panel represents the corresponding synthetic light curve 
in the 15-150 keV interval as a function of observer time t after 
To . The five coloured areas represent the integration times for the 
spectra in the lower panel. Except for the first interval they also 
correspond to the integration times of the optical measurements. 
The lower panel compares the five synthetic spectra to the multi 
wavelength data. The JHK fluxes were obtained during the last 
time interval where the burst was not detected in the optical 
ha.TiH. 
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Figure 9. Same as Fig. 8 with deceleration by a wind-like exter- 
nal medium (A* =0.1). 



erated electrons in internal shocks. However, we point out 
that, due to the long duration of this burst, it is highly 
probable that the deceleration of the outflow by the exter- 
nal medium will start before the end of internal shocks. 
Our analysis shows that a more realistic scenario taking 
into account this effect still agrees with the observations. 
The observed gamma-ray and optical flux would then be a 
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combination from internal and reverse shock contributions. 
Our modelling has been performed with an adopted redshift 
z = 0.3. Similar results would be obtained for a different z in 
the interval 0.05-0.85 (see section 3.3.1) with a suitable ad- 
justment of the model parameters (iJkin and micro-physics 
parameters) . 

The late near infrared data are more puzzling. The ob- 
served variability appears incompatible with a forward shock 
origin. On the other hand, the measured spectral slope is 

difficult to reproduce by the contribution of late internal 
shocks or the reverse shock, except if one assumes that the 
micro-physics parameters in the outflow can vary signifi- 
cantly during propagation. Finally, the early episodes in the 
gamma-ray light curve arc not compatible with most models 
of precursor activity (pliotosplioric emission, shock break- 
out, etc.) due to the very long delay before the main event. 
This would point out to a central engine having several dis- 
tinct episodes of relativistic ejection. 
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